Measurement of Inclusive Spin Structure Functions of the Deuteron 
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We report the results of a new measurement of spin structure functions of the deuteron in the 
region of moderate momentum transfer (Q 2 = 0.27 - 1.3 (GeV/c) 2 ) and final hadronic state mass in 
the nucleon resonance region (W = 1.08 - 2.0 GeV). We scattered a 2.5 GeV polarized continuous 
electron beam at Jefferson Lab off a dynamically polarized cryogenic solid state target ( 5 NDs) and 
detected the scattered electrons with the CEBAF Large Acceptance Spectrometer (CLAS). From 
our data, we extract the longitudinal double spin asymmetry An and the spin structure function 
gf. Our data are generally in reasonable agreement with existing data from SLAC where they 
overlap, and they represent a substantial improvement in statistical precision. We compare our 
results with expectations for resonance asymmetries and extrapolated deep inelastic scaling results. 
Finally, we evaluate the first moment of the structure function gf and study its approach to both 
the deep inelastic limit at large Q 2 and to the Gerasimov-Drell-Hearn sum rule at the real photon 
limit (Q 2 — > 0). We find that the first moment varies rapidly in the Q 2 range of our experiment and 
crosses zero at Q 2 between 0.5 and 0.8 (GeV/c) 2 , indicating the importance of the A resonance at 
these momentum transfers. 

PACS numbers: 13.60.Hb, 13.88.+e, 14.20.Dh 



I. INTRODUCTION 

The nucleon spin structure functions g p,n (x) and 
02 ( x ) an( i their moments have been extensively studied 
over the past two decades [llllllllllElillll[l(|. At 

large momentum transfer (Q 2 » 1 (GeV/c) 2 ) and final 
state mass (W > 2 GeV), these data can be successfully 
described via perturbative QCD (pQCD) up to next-to- 
leading order (NLO) and give us access to the helicity- 
weighted distribution functions Aq(x) and AG(x) of 
quarks and gluons in the nucleon [H], 0, 0, In 
this kinematic regime, one can relate the first moments 
= <?f ' (x)dx of the spin structure functions gi(x) 
(N = p or n) to the fraction of the nucleon spin carried 
by the quark hclicities and, via the famous Bjorken sum 
rule [l5lll6| . to the weak axial form factor gA- 

At lower momentum transfers, Q 2 1 (GeV/c) 2 , cor- 
rections proportional to powers of 1/ Q 2 develop due to 
higher twist and target mass effects [2> HE El m & d- 
dition to the logarithmic Q 2 dependence predicted by 
pQCD. As Q 2 decreases, an increasing part of the kine- 
matic range x = — 1 lies in the region of resonant final 
states (W < 2 GeV), which begin to dominate the spin 
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structure functions. They become less positive (or more 
negative in the case of the neutron), in particular in the 
region of the A resonance. Data in this region on struc- 
ture functions and on the (virtual) photon asymmetries 
Ai and Ai for the proton and the neutron, 

cr 1/2 -er 3 / 2 gi-gi/r 

Al = ; = — ? — ' 

C1/2 + C3/2 ^1 

A 2 = ° LT = ?lpl, (1) 
C1/2 + "3/2 V tF i 

can help us unravel the spin-isospin structure of reso- 
nance transition amplitudes and their interference with 
each other and with non-resonant terms. We can also test 
whether the observed duality between unpolarized deep 
inelastic and resonant structure functions ]2jJ,|2l|, |22| is 
realized for spin structure functions as well |23tl24j . Here, 
<7]/2 and (T3/2 are the (virtual) photon absorption cross 
sections for total (photon plus nucleon) helicity 1/2 and 
3/2 and <jlt is the longitudinal-transverse interference 
cross section, F\ is the unpolarized structure function, 
and r = v 2 IQ 2 with v = E — E' being the energy loss of 
the scattered electron. 

Due to the dominance of the resonances at low Q 2 , the 
integrals T p and T d m (TP + T n )/2 (which are positive 
in the scaling region of high Q 2 ) decrease rapidly and 
become negative as Q 2 approaches zero. In the limit 
Q 2 — > 0, the first moments for the proton and the neutron 
are constrained by the Gerasimov-Drell-Hearn Sum Rule 



3 



[HEl, which predicts that 




Here, a is the fine structure constant and M and kn 
are the mass and anomalous magnetic moment of the 
nucleon, respectively. Since the GDH sum rule is nega- 
tive, the integrals T\' d {Q 2 ) must have a negative slope at 
Q 2 = and then change rapidly at low Q 2 to meet the 
positive experimental results in the DIS region. 

So far, only phenomenological models fo r Tj (Q 2 ) co ver- 
ing the whole range of Q 2 exist pil^l^lMlMlsllsll. 
These models are constrained to reach the large-Q 2 
asymptotic value of the integral as measured by deep 
inelastic data and to approach zero at the photon point 
with a slope given by the Gerasimov-Drell-Hearn sum 
rule, Eq. J2J. The authors of Refs. jH,^ use a simple 
parametrization of the integral r 1+2 = J [gi(x) + gz{x)\dx 
to interpolate between these two points, and then sub- 
tract the integral over qg w hich is given by the Burkardt- 
Cottingham sum rule [34| . The approach taken in Refs. 
[30l l3ll uses a parametrization of existing resonance 
data and a vector meson dominance (VMD) inspired in- 
terpolation of the remaining integral strength at the two 
endpoints. 

For a complete picture of spin structure functions and 
their moments, one needs information on both the proton 
and the neutron. Since free neutron targets are imprac- 
tical, deuterium (as in the experiment described here) or 
3 He targets are used instead. An unambiguous extraction 
of neutron spin structure functions from nuclear ones is 
less straightforward in the resonance region than in the 
deep inelastic regime; however, the integrals are much 
less affected by uncertainties from Fermi motion, off-shell 
effects, and other nuclear corrections |35ll36ll37| . In par- 
ticular, studies j3^,|3jj show that the integral Tf for the 
deuteron from pion threshold on up is very close to the 
incoherent sum of proton and neutron integrals, once a 
correction for the deuteron D-state has been applied. 

So far, only very limited spin structure function data 
exist in the region of low to moderate Q 2 and W |40ll4l| . 
especially on the deuteron. A large program is underway 
at Jefferson Lab to map out the entire kinematic region 
Q 2 « 0.05 - 5 (GeV/c) 2 and W < 3 GeV. This pro- 
gram consists of measurements on 3 He (in Hall A) and 
on proton and deuteron targets with CLAS (the EG1 col- 
laboration in Hall B). First results from CLAS 01 an d 
Hall A |4^| have already been published. 

In the present paper, we present results on the 
deuteron from the first EG1 run in 1998, in which we 
measured double spin asymmetries A\\ = D{A\ + r\A-i) 
on deuterium with a beam energy of 2.5 GeV. (D and rj 
are kinematical factors, see Sec. II VI ) These data cover 
a range in Q 2 from 0.27 - 1.3 (GeV/c) 2 and final state 
mass in the resonance region (W = 1.08 - 2.0 GeV). The 
remaining data set from EG1 is presently under analysis 
and will increase both the kinematic coverage and the 



statistical precision of our data significantly. 

In the following, we give some details on the ex- 
periment (Sec. [Hj and its analysis (Sec. 1111(1 . We 
present our results on the deuteron spin asymmetry 
(Af +r)Af) (W, Q 2 ), the structure function gf(x,Q 2 ) 
and its first moment Tf (Q 2 ) (Sec. lIV|) . and conclude with 
a summary and outlook fSec. Iv|) . 

II. EXPERIMENTAL DETAILS 

The data described in this paper were collected during 
a three-month run in 1998, as part of the EG1 run group 
in Jefferson Lab's Hall B. A polarized electron beam with 
2.5 GeV beam energy was scattered off a deuterated am- 
monia ( 15 ND 3 ) target that was dynamically polarized 
along the beam direction. The average beam current 
of 2.5 nA corresponded to an instantaneous luminosity 
of 0.4 x 10 34 cm _2 s~ 1 . The beam polarization was mea- 
sured periodically with a M0ller polarimeter and the av- 
erage beam polarization was 72%. 

We used the CEBAF Large Acceptance Spectrometer 
(CLAS) to detect the scattered electrons. The CLAS de- 
tector |44| is built around six superconducting coils that 
produce a toroidal magnetic field. The orientation of the 
magnetic field can be chosen so that electrons are bent ei- 
ther toward (inbending) or away from the beam line (out- 
bending). The target was placed 55 cm upstream from 
its normal location in the center of CLAS to lower the an- 
gular threshold for electron detection and thus decrease 
the lower limit on the momentum transfer. Inbending 
electrons were detected down to a minimum polar angle 
of 14°. During this experiment the geometry of the tar- 
get excluded particle tracks with a polar angle between 
50° and 75°. The phi-acceptance is approximately 85%, 
limited mainly by the torus coils. 

The CLAS detector package consists of three layers 
of drift chambers for track reconstruction, one layer 
of scintillators for time-of-flight measurements, forward 
Cerenkov counters for electron-pion discrimination, and 
electromagnetic calorimeters to identify electrons and 
neutral particles. A coincidence between the Cerenkov 
and the calorimeter triggers the data acquisition. Elec- 
tron particle identification is accomplished using the 
Cerenkov detector and the distribution of energy de- 
posited in the calorimeter. The large acceptance of CLAS 
(« 1.5 sr for electrons) and its large kinematic coverage 
offset the limited luminosity that can typically be reached 
with polarized solid state targets (of order 10 35 cm _2 s _1 
at best), and allowed us to collect data for the entire W 
and Q 2 range simultaneously. 

The longitudinally polarized target was designed to fit 
within the 1 m central bore of the CLAS |45|. A pair 
of superconducting Hclmholtz coils provided a 5 T mag- 
netic field along the direction of the electron beam. The 
magnetic field was uniform to better than 1 x 10~ 4 in 
the center of the target over a length of 2 cm and a di- 
ameter of 2 cm. The ammonia crystals were contained 
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within a plastic cylindrical cell 1 cm in length and 1.5 cm 
in diameter. The cell was immersed in a liquid He bath 
maintained at approximately 1 K by a 4 He evaporation 
refrigerator. The cell was mounted on a target insert 
that also held a NH3 cell, as well as a 12 C and an empty 
cell. The latter cells were used to study the dilution of the 
measured asymmetries by events from unpolarized target 
constituents (see Sec. IIII 0*)! . The deuterons in the target 
were polarized using the Dynamic Nuclear Polarization 
(DNP) technique with 140 GHz microwaves. The 

polarization of the target was monitored on-line using 
the NMR technique. The NMR results were not used for 
our final analysis; instead, we extracted the product of 
beam and target polarization directly from our data, as 
described in Sec. IIII CI The beam was rastered on the 
ND3 target, although not over the full face of the target. 
The deuteron polarization suffered from this incomplete 
raster and from inadequate microwave power and ranged 
from approximately 10% to 25%. All data were taken 
with the target polarization along the beam direction, 
without reversal of the target polarization. The beam 
hclicity was reversed every second. 

During the 1998 run, we collected 300 million triggers 
for an integrated beam charge of about 0.4 mC. From 
this sample, 100 million electron events passed the cuts 
described in Sec. IIII Al These events covered a kinematic 
region from the quasielastic region (W « 0.94 GeV) to 
the edge of the deep inelastic region (W = 2 GeV) and 
for Q 2 = 0.27 - 1.3 (GeV/c) 2 . This kinematic coverage is 
shown in Fig. ^ together with the coverage of the second 
part of the EG1 experiment. 
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FIG. 1: Kinematic coverage of the data described in this 
paper ("EG1 2.5 GeV") together with the kinematic range of 
the second run of EG1 ("EG lb" at 1.6 GeV, 4.2 GeV, and 5.7 
GeV). The heavy solid lines indicate the elastic peak ( U W = 
0.9 GeV"), the location of the Su resonance ("W = 1.5 GeV") 
and the deep inelastic limit ( U W = 2 GeV" ) . Also shown are 
the kinematic lines for three representative values of x. 



III. DATA ANALYSIS 

The goal of our analysis was to first determine the 
double spin asymmetry 



*\\ = 



rU 



rll 



(3) 



for each kinematic bin and then to extract the physi- 
cal quantities of interest, the virtual photon asymmetries 
Af + r/A^ and the structure function gf, from the results. 
Here, er^-i stands for the differential electron scattering 
cross section with the target and electron spin pointing 
in opposite directions along the beam and correspond- 
ingly er^T f or parallel target and electron spin. 



A. Data selection 

For the present analysis, we selected data runs taken 
with a torus current of +2250 A (inbending electrons) 
and target polarization parallel to the beam direction. 
The data were taken with two slightly different beam en- 
ergies, 2.494 GeV and 2.565 GeV, due to a change of the 
accelerator configuration. We separated our sample into 
four different "run groups", two each with beam energy 
2.494 GeV and 2.565 GeV. Each run group corresponds 
to a contiguous set of runs with the same target material 
and approximately constant target parameters and run- 
ning conditions. Only runs with stable beam and detec- 
tor performance were included in our sample. The 2.565 
GeV groups also contained carbon target runs that were 
used to determine the dilution factor (see Sec. IIII 0]) . We 
analyzed events with scattering angles from about 14° to 
50° and scattered electron energies from 0.5 GeV to 2.5 
GeV. 

The data were sorted according to the helicity of the 
electron beam. During our run, the beam helicity fol- 
lowed a "pseudo-random" pattern of helicity pairs, where 
the first "bucket" (of 1 second length) of each pair was 
given random helicity and the second its complement. 
We matched the sequence of helicity bits for each event 
with the pattern sequence recorded in helicity scalers 
and discarded pairs for which the helicity assignment 
was inconsistent. We also discarded pairs with signifi- 
cantly different (by more than 10%) beam intensity in 
the two buckets (due to beam fluctuations or trips). The 
final data sample contained only matched pairs of buck- 
ets with stable running conditions. 

All events were accumulated in small bins of W 
(AW = 0.02 GeV) and Q 2 (AQ 2 /Q 2 w 20%), sepa- 
rately for both beam helicities. (The data on asymme- 
tries and gi, shown in Sec. II VI are weighted averages of 
several such bins). In addition, we also accumulated the 
integrated beam charge for each of the helicity buckets 
(corrected for deadtime) to normalize the helicity-sorted 
counts in each bin. We found that on average, there was 
a 0.3% difference between the integrated charge for the 
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two opposite helicities, possibly stemming from the sen- 
sitivity of the photocathodc in the polarized source to 
small remaining linear polarization components or beam 
motions of the photo-ionization laser beam. Our nor- 
malization method removed the effect of this asymmetry, 
and it was further suppressed by reversing the relative 
sign between the helicity at the cathode and in the ex- 
perimental Hall (through spin-precession in the injector 
and the accelerator). 



B. Electron cuts 

We selected electron events by first requiring a negative 
track with matching signals in the time-of-flight (ToF) 
scintillators, the Cerenkov counters (CC), and the elec- 
tromagnetic calorimeter (EC). In the presence of several 
such tracks, the track with the shortest flight time was 
selected as the electron candidate. Some additional cuts 
on the track vertex along the beam line removed events 
from the entrance and exit windows of the polarized tar- 
get chamber, as well as badly reconstructed tracks. 

We used information from the CC and the EC to fur- 
ther separate electrons from negative pions. We required 
a signal in the CC that exceeded 50% of the average 
signal for a single photo-electron. Furthermore, we re- 
quired that the energy measured in the EC exceeded 20% 
of the candidate electron momentum (the average sam- 
pling fraction of the EC was 27%). A typical example for 
the ratio of sampled EC energy over momentum is shown 
in Fig. |2 The open histogram shows events that passed 
all other electron cuts (including the CC cut). 

We also collected a sample of ir~ events with no signal 
above threshold in the CC. As shown by the shaded area 
in Fig. |21 the E/p spectrum associated with n~ events 
is strikingly different from the electron spectrum. Un- 
der the conservative assumption that all events below a 
E/p ratio of 0.15 came from pions, we cross-normalized 
the two spectra below that point and estimated the re- 
maining pion contamination of our electron sample by 
the ratio of the two integrated spectra above our cut of 
E/p > 0.2. For all kinematics studied, this remaining 
contamination turned out to be less than 1%. 

The reconstructed momenta of the scattered electrons 
were corrected for effects from unknown torus field distor- 
tions and slight drift chamber misalignments. We used 
NH3 runs taken interleaved with the ND3 ones to deter- 
mine the correction factor by optimizing the position and 
width of the elastic peak (W — 0.938 GeV) for all scat- 
tering angles 9 and <f>. The resulting corrections were of 
the order 0.1% on average. 



C. Dilution and polarization 

The double spin asymmetry An can be extracted from 
the count rate asymmetry (normalized by the integrated 




E/p 



FIG. 2: Spectra of the ratio of measured energy, E (in GeV), 
in the electromagnetic calorimeter over the track momen- 
tum, p (in GeV/c), for electrons (open histogram) and pions 
(shaded area). The vertical scale is arbitrary. Both spectra 
have been cross-normalized at low E/p. Events above the 
indicated threshold are identified as electrons. 

beam charge) after accounting for the dilution from unpo- 
larized target constituents and the beam (Pi,) and target 
(P t ) polarization: 

A meas = 1 N+/Q+-N~/Q- 

II DF-P b -P t N+/Q+ + N-/Q- 1 1 ' 

where N + ' ~ are the counts and Q +: ~ are the integrated 
beam charge for positive and negative helicity. 

We determined the dilution factor DF in Eq. (J2J by 
approximating the contribution to the count rates from 
unpolarized target constituents (target foils, LHe coolant 
and 15 N in ammonia) with the spectra taken on the car- 
bon target. Some components of these two targets were 
the same (e.g., the LHe coolant and foils were present 
for the carbon target as well) , and carbon, nitrogen, and 
even He have similar binding energies per nucleon and 
Fermi momenta, suggesting that their inclusive electron 
scattering spectra are similar after correcting for the total 
number of target nucleons. (This assumption has since 
been verified to better than 3% with dedicated runs on a 
pure 15 N target during the second part of EG1). 

To account for the different number of nucleons in each 
target and different overall target thicknesses, we cross- 
normalized the carbon target spectra to the ammonia 
target spectra. We determined a normalization constant 
A such that the two spectra had the same number of 
counts below a cut-off missing mass W cut , well below 
the quasi-elastic peak. The cut-off ranged from W cu t — 
0.835 GeV at Q 2 = 0.3 (GeV/c) 2 to W cut = 0.5 GeV at 
Q 2 = 1.2 (GeV/c) 2 , and was chosen so that the deuteron 
contribution was negligible, according to a Monte Carlo 
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simulation of the deuteron wave function. 



& 

Z5 

o 
O 

"D 

CD 
N 



16000 
14000 
12000- 
10000 
8000 
6000- 
4000 
2000 




i ND3 

• Carbon 

* Deuteron 

— Deuteron (quasi-elastic) 




0.75 




1.25 



1.5 



1.75 



W [GeV] 

FIG. 3: Spectra of counts vs. final state mass W from the 
polarized ND3 target (solid triangles) and carbon target (cir- 
cles) runs for the range Q 2 = 0.5±0.1 (GeV/c) 2 . The spectra 
have been cross-normalized at low W. The deuteron spectrum 
(open triangles) is the difference between these two spectra. 
In the quasi-elastic peak region, it agrees well with a simu- 
lation using the Paris wave function for the deuteron (solid 
line) . 



The dilution factor can then be written as 
N Nr > 3 -A-N c 



DF 



N. 



ND 3 



(5) 



where the numerator is the count rate due to deuterium 
alone. The results of this method for an intermediate Q 2 
bin are shown in Fig. [21 The normalized carbon spectrum 
(circles) has been subtracted from the ammonia spectrum 
(solid triangles) to yield the deuteron spectrum (open 
triangles). The line indicates the result of our Monte 
Carlo simulation of the deuteron spectrum alone, which 
is based on quasi-elastic scattering (PWfA) and the Paris 
wave function |48j for the deuteron. The dilution factor 
for our experiment was around DF « 0.2. 

The second ingredient needed in Eq. is the prod- 
uct of beam and target polarization. We measured both 
the beam polarization (with a M0ller polarimeter) and 
the target polarization (using NMR) individually during 
the run. However, due to the small amount of target 
material and its inhomogeneous exposure to the electron 
beam, the NMR results were not very precise and reli- 
able. Instead, we determined directly the product Pb ■ Pt 
by extracting it from the measured asymmetry in the 
quasi-elastic region. For this purpose, we used inclusive 
quasi-elastic events d(e, e') in the range 0.85 GeV< W < 
1.0 GeV. 

The asymmetry A| | for elastic scattering from protons 
and neutrons can be calculated from known nucleon form 
factors with very little systematic uncertainty (less than 



1-2% in our kinematic region). We used our simulation 
of the deuteron wave function to calculate the expected 
asymmetry for inclusive quasi-elastic scattering within 
our kinematic cuts, which differed only slightly from the 
cross section-weighted average of the proton and neutron 
asymmetries. We used the dilution factor determined by 
the method described above to extract the product Pb-Pt- 
Due to the large kinematic coverage of CLAS, data 
on the quasi-elastic asymmetries were collected contin- 
uously and simultaneously with the inelastic asymmetry 
data. The extracted average polarization product Pb ■ Pt 
for each of the four run groups is therefore a faithful rep- 
resentation of the running conditions for that group, with 
minimal systematic uncertainties. Our results are shown 
in Fig. where we divided the product Pb ■ Pt by the 
measured beam polarization to extract the target polar- 
ization. The results for each of the individual run groups 
have statistical errors on the order of 13%, which were 
included in the total statistical error of the asymmetries 
from each run group. The final results for the inelastic 
asymmetries are statistically weighted averages from the 
four run groups, with a contribution to their statistical 
errors from the polarization product of about 6.7% of 
their values. 




Run Group 

FIG. 4: Average target polarization for each of the four run 
groups determined by dividing the values of the product Pi, • 
Pt extracted from the quasi-elastic asymmetry by the beam 
polarization measured with a M0ller polarimeter. The target 
polarization decreases over time due to beam exposure. 



D. Other backgrounds 

After dividing out the dilution factor and the beam 
and target polarizations in Eq. Q), we corrected the ex- 
tracted asymmetry for additional background contribu- 
tions. These include contamination of the scattered elec- 
tron sample by negative pions and pair-produced elec- 
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trons, as well as contributions from polarized target con- 
stituents other than deuterium. 

We already discussed the contribution from pions 
misidentified as electrons, which was less than 1% in all 
cases. A more important contribution comes from elec- 
trons that are decay products of neutral pions (either 
through the Dalitz decay ir° — ► je + e~ or pair conver- 
sion of decay photons). The rate of electrons from these 
decays was estimated using the Wiser fit 0] for pion 
photoproduction and tested against the Monte Carlo 
code "PYTHIA" . We also measured directly the rate of 
positron production in each kinematic bin (again making 
use of the large acceptance of CLAS for both positively 
and negatively charged particles). This rate should be 
equal to that of electrons from charge-symmetric decays 
and was found to agree well with the Wiser fit. The 
asymmetry for positrons was found to be consistent with 
zero and in any case no larger than the asymmetry for 
electron scattering events. We used a parametrization of 
our results to estimate the fraction of detected electrons 
coming from these decays. This fraction was typically 
1% for most of the kinematic region, but increased up 
to 20% at the highest W values. We corrected our data 
for this background by applying a further dilution factor 
to our asymmetries. Since we could not exclude a small 
non-zero asymmetry for these events, we assumed a sys- 
tematic uncertainty equal to the size of this correction. 

The nitrogen in our dynamically polarized ammonia 
target carries a small residual polarization, which leads 
to a partially polarized bound proton in 15 N. Possible 
additional polarized target species include isotopic impu- 
rities of 1 N and X H. Extensive experience with similar 
targets at SLAC 0] shows that the corresponding correc- 
tions to the asymmetry are at most a few percent. We 
included the uncertainty due to these contributions in 
our systematic error. 

Another potential contribution to the measured asym- 
metry comes from parity-violating electron scattering off 
all target constituents. However, at the low momentum 
transfers of our experiment, the expected asymmetry is 
less than 10" and can be treated as another (small) 

systematic uncertainty. 



E. Radiative corrections and models 

The final step in the extraction of the desired "Born 
asymmetry" A\\ requires correcting the measured asym- 
metry for higher-order electromagnetic processes (in- 
ternal radiative corrections) and electron energy loss 
through bremsstrahhmg in the target before or after the 
scattering (external radiative corrections). These radia- 
tive corrections were applied separately to the numerator 
and the denominator of Eq. @ , which yields an additive 
{■Arc) and a multiplicative (Frc) correction term: 

A^A^/Frc + Arc. (6) 



Here, the factor 1/Frc represents the increase of the 
denominator in Eq. due to the radiative elastic and 
quasi-elastic tails that act like an additional dilution of 
the inelastic events. Correspondingly, the statistical error 
of the final result was scaled up by 1 /Frc as well. 

Both components (Frc, Arc) were determined by 
running the code "RCSLACPOL" developed at SLAC |jf . 
This code uses parametrizations of all relevant input 
quantities (structure functions and form factors), as well 
as a model of our target, to calculate both fully radiated 
and Born cross sections and asymmetries. It is based on 
the approach developed by Kuchto and Shumeiko pHl for 
the internal corrections and by Tsai |H2] for the external 
corrections, including the radiative depolarization of the 
beam due to external bremsstrahlung. 

We used parametrizations of the world data on polar- 
ized and unpolarized structure functions and elastic form 
factors as input for the radiative correction code and to 
extract physics quantities of interest from the measured 
asymmetries. These parametrizations are described in 
and are based on fits to unpolarized structure function 
data from NMC [Hjj and SLAC [HI IM IHl l57|and polar- 
ized structure function data from SL AC0,|l|i g HJ , 
CERN 0,013, and HERMES MJI3- The nucleon form 
factors were taken from Ref. p3 with updated values 
for the ratio Ge p /Gmp from the recent Jefferson Lab ex- 
periment |59|. For the asymmetries A\ and Ai in the 
resonance region, we used parametrizations of resonance 
transition amplitudes from Ref. .3Q| (in the form of a 
computer code named "AO") and Ref-pJ ("MAID") 
together with a fit of the SLAC data |4l|. We also in- 
cluded our own preliminary asymmetry data in these fits. 
All fits were varied within reasonable errors or replaced 
with alternative existing fits to study the systematic de- 
pendence of our final results on these parametrizations. 



F. Systematic errors 

The total systematic error on our data ranges from 
25% to 50% of the statistical error for the asymmetries 
and from 35% to 50% of the statistical error for the 
structure function gf. The leading contributions to these 
systematic uncertainties come from radiative corrections 
(40-50% of the total systematic error on average) , uncer- 
tainties in the unpolarized structure functions needed to 
extract final physics results (also 40-50% of the total), 
and the dilution factor (about 40%). We also considered 
the effect of finite resolution and errors in the measured 
kinematic variables (about 10% of the total). At higher 
Q 2 and especially higher W , pair-symmetric decay elec- 
trons also contributed significantly to the overall system- 
atic uncertainty (15-20% averaged over all kinematic bins 
and most of the systematic error at the kinematic limit). 
Finally, for the extraction of the spin structure function 
gf and its integrals, some model assumption about the 
virtual photon asymmetry A2 is needed (see Sec. lIVf) and 
leads to a further systematic error (up to 50%). 
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We accounted for each of these systematic errors by 
changing a relevant input parameter or model, and then 
repeating the entire analysis up to the final results, in- 
cluding the integrals of gf over the measured region. We 
took the error as the deviation of the alternative results 
from the standard analysis. We added all uncorrelated 
systematic errors in quadrature. The final systematic er- 
rors are shown in the data tables in Sec. II VI 

For the radiative correction errors, we varied all in- 
put models and parametrizations for the radiative code, 
including polarized and unpolarized structure functions, 
form factors, and the target model, within realistic lim- 
its. We also checked the accuracy of the peaking approx- 
imation by comparing the results with those from a full 
integration without approximations. 

Similarly, we varied the models for the unpolarized 
structure functions F? and R = — , which entered the 

extraction of gf and the asymmetry Af + r]Af; from our 
data (see Sec . IIVII . We used different fits of the world 
data [54L 1571 l61| . and studied their effect on the final 
physics results. In the case of the polarized structure 
function gf and its integrals, we also varied the model 
for the asymmetry Af; from Af = to the prediction by 
the "MAID" code and a simple parameterization based 
on the twist-2 result by Wandzura and Wilczek [fj^] that 
describes the SLAC data @ well. 

For the error introduced by the uncertainty in the di- 
lution factor, we varied the cross-normalization between 
the carbon and ammonia target data by an amount of 
6%, consistent with the variations observed for different 
W and Q 2 ranges and possible differences in the 12 C and 
15 N spectra. This yields an average variation of the di- 
lution factor Eq. JSJ of 25%, making this error a safe 
upper bound for all systematic errors that are directly 
proportional to the measured asymmetry. 

The CLAS momentum resolution and reconstruction 
effects were studied by moving all data points by 0.02 
GeV in W and recalculating the final results. The effect 
of this variation on the integrals of gf also gave an upper 
limit to systematic errors due to the integration method, 
which consisted in a simple sum of all bins of size AW — 
0.02 GeV, multiplied by the bin width in x. 

Other systematic errors were either negligible or have 
already been described in the previous section. We note 
that we do not have a significant systematic error from 
the beam and target polarization product, since they 
were directly determined from our data (with minimal 
theoretical uncertainty). In particular, the theoretical 
asymmetry A^ as is only weakly dependent (±1%) on the 
elastic form factor ratio Ge/Gm for the proton. How- 
ever, the statistical error of this method is not negligible 
and was included in the total statistical error of the final 
results. 



1.0 




1.0 1.2 1.4 1.6 1.8 2.0 2.2 

W [GeV] 



FIG. 5: Af + r]Ai versus W for Q 2 = 0.39 - 0.65 (GeV/c) 2 . 
Our data points are indicated as triangles with statistical er- 
rors only. The size of the systematic error is indicated by the 
shaded band at the bottom of the graph. Previous data from 
SLAC E143 |fj are shown as open circles with statistical and 
systematic errors combined. The positions of several promi- 
nent resonances are indicated by the labeled arrows. The solid 
line is our model parametrization of the world data (without 
nuclear corrections like Fermi motion and off-shell effects) 
and the dashed line is the resonant contribution to Af alone 
(from the code "AO"). 



IV. RESULTS 

A. Virtual photon asymmetries 

We extracted a combination of the virtual photon 
asymmetries, Af + rjAf, from our data on An using a 
parametrization |57| of the structure function R, via the 
relationship 

An=D(A 1 +r l A 2 ), (7) 

where the virtual photon depolarization factor is given 
by D = (1 - eE'/E)/(l + eR) and rj = c*JCp/(E - eE') 
(e is the virtual photon polarization parameter, E is the 
beam energy and E' is the scattered electron energy). 

The extracted photon asymmetries (Af + rjAf) (W, Q 2 ) 
for three different Q 2 bins are listed in Tables HI IIIII to- 
gether with their statistical and full systematic errors. 
We show the results for our intermediate Q 2 bin in Fig. EI 
together with previous data from SLAC || and some 
model calculations. A comparison of the three different 
Q 2 bins can be found in Fig. 

Since we didn't measure the asymmetry with the target 
polarization perpendicular to the electron beam ( A± ) , we 
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TABLE I: The measured virtual photon asymmetry A 1 
of the deuteron for Q 2 = 0.27 - 0.39 (GeV/c) 2 . 



VV [uc V J 


A d -A- n A d 


OLdL. HiLUJi 


Syst. Error 


1 in 


A. QflO 


U.ooU 


n ori7 


i on 


n 07Q 


U.zUo 


n aai 
U.UD1 


1 OS 
l.Zo 


-U.^UO 


u.ioy 


U.Uol 


i oa 
1.36 


A OOO 


u.iyi 


0.069 


1.44 


-0.124 


0.161 


0.028 


1.52 


-0.077 


0.131 


0.017 


1.60 


-0.036 


0.119 


0.015 


1.68 


0.140 


0.102 


0.023 


1.76 


0.063 


0.101 


0.011 


1.84 


0.055 


0.086 


0.017 


1.92 


-0.254 


0.080 


0.028 


2.00 


-0.084 


0.072 


0.009 



TABLE II: The measured virtual photon asymmetry A 1 

0.39 - 0.65 (GeV/c) 2 . 



rjAi of the deuteron for Q 2 



W [GeV] 


Ai + nAi 


Stat. Error 


Syst. Error 


1.12 


-0.327 


0.267 


0.191 


1.20 


-0.411 


0.109 


0.081 


1.28 


-0.316 


0.090 


0.061 


1.36 


-0.070 


0.101 


0.062 


1.44 


0.086 


0.085 


0.022 


1.52 


0.144 


0.068 


0.025 


1.60 


0.147 


0.063 


0.024 


1.68 


0.061 


0.054 


0.015 


1.76 


0.006 


0.053 


0.011 


1.84 


0.024 


0.050 


0.013 


1.92 


-0.045 


0.047 


0.013 



cannot directly extract the asymmetry Af or Af,- The in- 
terference term A 2 is limited by \A 2 \ < \J R(A± + l)/2, 
where the value of R is around 0.1 - 0.3 at Q 2 = 0.5 
(GeV/c) 2 H3 and the typical size of r\ for our experi- 
ment ranges from 0.1 at W = 2 GeV to 1.2 right at pion 
threshold (W = 1.08 GeV). Correspondingly, the asym- 
metry A% could contribute as much as 0.07 (high W) to 
0.15 (at threshold) to the asymmetries shown in Figs.[S] 
and|S| However, according to our paramctrization, this 



TABLE III: The measured virtual photon asymmetry Af 
r)Ai of the deuteron for Q 2 = 0.65 - 1.3 (GeV/c) 2 . 



W [GeV] 


Af + n A d 2 


Stat. Error 


Syst. Error 


1.12 


-0.529 


0.223 


0.125 


1.20 


-0.299 


0.101 


0.038 


1.28 


-0.106 


0.083 


0.025 


1.36 


-0.005 


0.091 


0.046 


1.44 


0.139 


0.078 


0.017 


1.52 


0.340 


0.067 


0.035 


1.60 


0.307 


0.061 


0.038 


1.68 


0.195 


0.054 


0.027 


1.76 


0.184 


0.056 


0.033 



< 




FIG. 6: Our data for three different bins in Q 2 , together with 
statistical errors. (Systematic errors are highly correlated be- 
tween different Q 2 bins and should have only minor effects on 
the observed Q 2 -dependence). The long-dashed line shows 
our model parametrization of Af+r^A^ for Q 2 — 1.0 (GeV/c) 2 
and the short-dashed line shows our model for Q 2 = 0.34 
(GeV/c) 2 . 



contribution should be more typically of order 0.02. 

With this caveat, one can conclude that the data shown 
in Fig.[S]exhibit the expected behavior for the asymmetry 
Af. In the region of the A(1232) resonance, the asym- 
metry is strongly negative, and fully compatible with 
the expectation Af = —0.5 for the resonance contri- 
bution alone. Beyond W = 1.4 GeV, the asymmetry 
becomes positive, indicating that helicity-1/2 transition 
amplitudes begin to dominate even at this rather low 
Q 2 . However, even in the region of the Sn resonance the 
asymmetry is markedly smaller (around 0.15) than for 
the proton (around 0.5, see Ref. yj), indicating that for 
the neutron alone the helicity-3/2 amplitude may still 
be larger. Figure also shows the predicted full asym- 
metry from our parametrization and a prediction for the 
resonance contributions to Af alone. The latter is based 
on the code "AO" [3(j , which uses a fit of exclusive pion 
electro- and photoproduction data to parametrize reso- 
nant and Born pion production amplitudes. Apparently, 
the contribution from the resonances alone already de- 
scribes the data well in the region of low to intermediate 
W , while non-resonant contributions (and maybe a siz- 
able asymmetry Af) are needed at high W. In general, 
our data agree fairly well with model predictions and the 
existing SLAC data. However, they have significantly 
smaller statistical errors and better resolution in W, as 
well as coverage down to lower Q 2 than the SLAC data. 
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TABLE IV: The spin structure function gf of the deuteron 
for Q 2 = 0.27 - 0.39 (GeV/c) 2 . 



W [GeV] 


gf 


Stat. Error 


Syst. Error 


1.12 


0.033 


0.058 


0.042 


1.20 


-0.115 


0.080 


0.029 


1.28 


-0.172 


0.074 


0.033 


1.36 


-0.080 


0.067 


0.022 


1.44 


-0.059 


0.067 


0.015 


1.52 


-0.040 


0.078 


0.010 


1.60 


-0.022 


0.073 


0.013 


1.68 


0.112 


0.074 


0.020 


1.76 


0.058 


0.075 


0.011 


1.84 


0.048 


0.065 


0.015 


1.92 


-0.202 


0.066 


0.028 


2.00 


-0.070 


0.066 


0.012 



TABLE V: The spin structure function gf of the deuteron 
for Q 2 = 0.39 - 0.65 (GeV/c) 2 . 



T XT „t ri 


„d 

9i 


Stat. Error 


Syst. Error 


1.12 


-0.025 


0.018 


0.004 


1.20 


-0.088 


0.025 


0.017 


1.28 


-0.083 


0.024 


0.011 


1.36 


-0.008 


0.023 


0.014 


1.44 


0.024 


0.024 


0.007 


1.52 


0.075 


0.029 


0.012 


1.60 


0.080 


0.028 


0.014 


1.68 


0.046 


0.030 


0.012 


1.76 


0.016 


0.031 


0.011 


1.84 


0.028 


0.030 


0.013 


1.92 


-0.016 


0.031 


0.013 



A comparison of our results for different Q 2 (see Fig. El) 
shows a general trend toward more positive asymmetries 
for higher Q 2 , especially in the region of the Sn and Dn 
resonances. This is in agreement with the expected tran- 
sition from helicity-3/2 dominance at low Q 2 (and espe- 
cially at the photon point, where it yields the negative 
value for the GDH sum rule), and helicity-1/2 dominance 
at higher Q 2 . In the limit of very large Q 2 , the asym- 
metry Af in the resonance region should become close to 
1, as predicted by pQCD as well as hyperfine-improved 
quark models and duality arguments. A similar behavior 
is observed for the proton asymmetries 



TABLE VI: The spin structure function gf of the deuteron 
for Q 2 = 0.65 - 1.3 (GeV/c) 2 . 



W [GeV] 


d 

9i 


Stat. Error 


Syst. Error 


1.12 


-0.022 


0.008 


0.004 


1.20 


-0.029 


0.011 


0.003 


1.28 


-0.010 


0.011 


0.004 


1.36 


0.003 


0.011 


0.007 


1.44 


0.024 


0.013 


0.005 


1.52 


0.086 


0.016 


0.011 


1.60 


0.089 


0.016 


0.013 


1.68 


0.072 


0.018 


0.013 


1.76 


0.082 


0.021 


0.017 



B. Spin structure function gf 

The spin structure function gf(W, Q 2 ) was calculated 
from the photon asymmetry (Af + r)A^){W, Q 2 ) for each 
bin using 

9t (W, Q 2 ) = j^(Ai + j-A*) Ff(W, Q 2 ) 

Here, Ff w (Ff + F[ l )/2 represents the unpolarized 
structure function of the deuteron (per nuclcon) and 
r = v 2 IQ 2 . Because of the partial cancellation of the 

two terms in ( — r^j,gf is less sensitive to the asym- 
metry A-i . We list our results for gf with their statistical 
and full systematic errors (including the uncertainty due 
to A 2 ) in Tables HH- EH 

In Fig. [7J we show our results for all three values 
of Q 2 , plotted against the Nachtmann scaling variable 
£ = Q 2 /M(v+q). This variable corresponds to Bjorken x 
at high Q 2 while it takes target nucleon mass corrections 
into account and therefore reduces "kinematical higher 
twist" scaling-violating effects at lower Q 2 . Together 
with our data, we also show as reference the prediction 
for gf(£,Q 2 — 5 (GeV/c) 2 ) from our model. The assump- 
tion of local quark-hadron duality predicts that structure 



functions like F± and g\ should, on average, approach a 
universal scaling curve if plotted versus the variable £, 
even in the resonance region. This is confirmed down to 
rather low Q 2 in the case of the unpolarized structure 
function F% 0, [2i| • Apparently, local duality does not 
work as well for the polarized structure function gf at 
high values of £ where the asymmetry is dominated by 
the A resonance and therefore is negative. Overall, the 
approach to the "asymptotic value" for Q 2 — 5 (GeV/c) 2 
seems to be relatively slow; only our highest Q 2 bin shows 
fairly good agreement beyond the region of the A reso- 
nance. 



C. Integrals 

We calculated the integrals Tf(Q 2 ) = J gf(x,Q 2 )dx 
for our results on gf(x,Q 2 ) over the (ordinary) Bjorken 
variable x for four different Q 2 bins, beginning at quasi- 
free pion production threshold (W = 1.08 GeV) up to 
the kinematic limit of our data. (The first two Q 2 bins 
are the same as shown in Tables IIVI and while we 
split the last bin into two halves, from Q 2 = 0.65 to 0.92 
(GeV/c) 2 and from Q 2 = 0.92 to 1.3 (GeV/c) 2 .) We 
expect that these integrals are close to an incoherent av- 
erage over the individual nucleons (proton and neutron) 
in deuterium, reduced by the D-state correction factor 
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FIG. 7: The spin structure function gf for the deuteron at 3 
different values of Q 2 , plotted against the Nachtmann variable 
£ together with an extrapolation of a fit to the deep inelastic 
data at Q 2 = 5 (GeV/c) 2 . Following standard conventions, all 
values are normalized to the number of nucleons in deuterium. 
The error bars are statistical only, while the shaded bands 
indicate systematic error bars for the three data sets. 



TABLE VII: The first moments of the spin structure function 
gf of the deuteron. Following standard convention, the inte- 
gral is normalized to the number of nucleons in deuterium. 
Q 2 is in (GeV/c) 2 and W max in GeV. 



Q 2 


Wmax 


Meas. Ti 


Stat. Err. Sys. Err. 


Full Ti 


Sys. Err. 


0.34 


2.00 


-0.027 


0.012 


0.005 


-0.034 


0.008 


0.53 


2.00 


-0.008 


0.004 


0.002 


-0.013 


0.007 


0.79 


1.96 


0.008 


0.003 


0.003 


0.009 


0.008 


1.10 


1.80 


0.007 


0.003 


0.002 


0.016 


0.009 



(1 — 1.5Pd), where Pd ~ 0.05 is the deuteron D-state 
probability. The results are shown in the third column 
of Table IV11I and the upper kinematic limits for W are 
listed in the second column. These upper W bounds cor- 
respond to lower limits of x — (0.1,0.15,0.21,0.32) for 
the four Q 2 bins, respectively. 

We use our model to estimate the contribution to 
the integral below these limits and show the resulting 
"full" integrals and their systematic errors in the last two 
columns of Table IVIII These systematic errors include a 
contribution from the uncertainty of this extrapolation to 
x = 0. To estimate this uncertainty, we studied the vari- 
ation of the low-x contribution according to different fits 
to the world data; also, since there are few high-precision 
data below x = 0.03, we added a systematic error equal 
to the value of the integral from x — up to 0.03. Due 
to the large theoretical uncertainty about the shape of 



the spin structure functions at very low x and the ab- 
sence of high-precision data in this region, the error on 
this extrapolation may be even larger than indicated by 
our systematic error estimate (see below). 
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FIG. 8: The first moment of the spin structure function g 1 
of the deuteron (per nucleon). See explanations in text. 

Our results for the first moment Tf (Q 2 ) of the spin 
structure function gf are shown in Fig. |S| The solid line 
at higher Q 2 is a fit to the world's data in the DIS region 
including QCD corrections up to second power in the 
strong coupling constant. The dotted line indicates the 
slope for the integral at Q 2 = predicted by the GDH 
sum rule (we use the incoherent sum of the results for 
the proton and for the neutron, normalized to two). The 
short-dashed line is the result from the code "AO" [3(j 
for the contribution from the nucleon resonances only. 
The long-dashed line by Burkert and Ioffe 0, is 
the "AO" result plus a term that depends smoothly on 
Q 2 and interpolates between the part that is missing at 
Q 2 = to saturate the GDH sum rule and the full value 
of Ti in the high Q 2 limit. Fig. |S| also shows the pre- 
diction from the model by Soffer and Teryaev [2^, [2j| 
(dot-dashed line). They use an interpolation of the in- 
tegral over the structure function g T — g x + gi, which 
converges to Ti at high Q 2 and remains positive down 
to the photon point where its slope is given by a com- 
bination of the nucleon charge and anomalous magnetic 
moment. They subtract the contribution from the inte- 
gral over gi (which is related to nucleon form factors via 
the Burkardt-Cottingham sum rule) to obtain the inte- 
gral Ti alone. The same authors have recently published 
a new parametrization of the proton-neutron difference 
integral for all Q 2 (6^ which might change the curve for 
the deuteron shown here. The solid triangles are based 
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on EG1 data alone and the open triangles include the 
estimated contribution to the integral from beyond our 
kinematic limits. The inner error bars are statistical and 
the outer error bars represent the systematic errors added 
in quadrature. They include the uncertainty on the es- 
timated low- a: contribution for the full integrals (open 
triangles). 

The first conclusion one can draw from Fig. [S] is that 
the integral over our measured region (essentially the res- 
onance region) is in rather good agreement with the pre- 
diction of the "AO" parametrization for resonance con- 
tributions only. The data follow the predicted trend from 
negative values at small Q 2 , where the A resonance con- 
tributes most of the integral and most other resonances 
are also dominated by the helicity-3/2 transition ampli- 
tude, to positive values at higher Q 2 , where the helicity- 
1/2 amplitude begins to take over and the importance of 
the A is diminished. Since we did not include Born terms 
or other non-resonant terms in the curve labeled "AO" , 
one can conclude that these terms must contribute rela- 
tively little to the integral over the resonance region in 
the case of the deutcron. This may be due to a partial 
cancellation between the asymmetry of the proton (which 
is likely positive for these terms) and that of the neutron. 

Extrapolating the integral down to x — seems to 
change the results only moderately (in the negative di- 
rection at low Q 2 and towards more positive values at 
higher Q 2 ). This can be understood again as a cancella- 
tion between a strongly negative-going trend of the struc- 
ture function g™(x) as x goes to zero and a more positive 
trend for of (xb according to existing DIS data and NLO 
analyses |l4ll64|. However, at present our understanding 
of the behavior of spin structure functions at very low 
x is still incomplete, making this extrapolation rather 
uncertain (as it is in the DIS region). Therefore, the er- 
ror bars on our open triangles may still underestimate 
that uncertainty. The emergence of new information on 
the low-a; behavior of spin structure functions over the 
past five years is responsible for most of the apparent 
disagreement between our quoted results and those from 
the E143 experiment at SLAC. The integrals over the 
resonance region alone agree fairly well with the SLAC 
data (to within 1.1 standard deviations); however, the ex- 
trapolation beyond W = 2 GeV is much more negative 
for the parametrization used in the present analysis and 
would move the SLAC data points down by about 0.008 
and 0.015 at Q 2 = 0.5 and 1.1 (GeV/c) 2 , respectively 
With this proviso, our data are (marginally) consistent 
with the SLAC data, but have much improved statistical 
errors and cover lower Q 2 . 

Our data lie somewhat below both phenomenological 
predictions for the full integral shown in Fig. 03 suggest- 
ing a slower transition from the negative values near the 
photon point to the positive asymptotic value at high Q 2 . 
The zero-crossing appears to occur somewhere between 
Q 2 = 0.5 (GeV/c) 2 and Q 2 = 0.8 (GeV/c) 2 , significantly 
later than in the case of the proton |£g. However, the 
systematic errors are highly correlated point-to-point so 



that the deviation from the predictions b y B urkert and 
Ioffe HI i3 and by Soffcr and Teryaev jH is not 
highly significant. 



V. SUMMARY AND OUTLOOK 

In this paper, we report first results on deuterium for 
inclusive spin structure functions in the nuclcon reso- 
nance region from the EG1 program at Jefferson Lab. 
These data significantly expand the kinematic coverage 
and statistical precision beyond the only previous data 
from SLAC [4lJ . We find generally reasonable agreement 
between these two data sets and various model predic- 
tions and parametrizations. In particular, the impor- 
tance and the negative asymmetry of the A resonance is 
confirmed, as is the general trend to more positive asym- 
metries at higher Q 2 and W. 

The spin structure function gf is less positive than for 
the proton case indicating that the neutron contri- 
bution is mostly negative in our kinematic region. While 
9i (£j Q 2 ) seems to be approaching the DIS scaling curve 
for large W and Q 2 , there are significant deviations from 
"local duality" , again mostly due to the A resonance. 

The integral over gf follows the expected trend in gen- 
eral, rising towards the DIS limit at the highest measured 
Q 2 while dropping rapidly below towards our lowest Q 2 
point. Clearly, neither the kinematic reach (in W and 
Q 2 ) nor the statistical precision of the present data set 
allow a definite statement about the validity of (or the ap- 
proach towards) the GDH sum rule limit. However, our 
data constrain the general trend required of any theory 
that aims to describe the spin structure of the nucleon 
over the full range of length scales, from the real photon 
point to the scaling limit. 

Spin structure function data on the deuteron, together 
with the corresponding proton results, should in principle 
allow us to separate the different isospin contributions to 
the resonant and non-resonant asymmetries. However, 
the first run of EG1 analyzed here did not yield enough 
statistical precision to make a direct separation of pro- 
ton and neutron contributions to the deuteron asymme- 
try feasible. However, we plan to submit results on the 
integral Ti for the neutron and the proton-neutron dif- 
ference, extracted from our data on the proton and the 
deuteron, in a separate paper. In the meantime, the com- 
plete EG1 data set has been collected in a second run, 
which will yield a nearly tenfold improvement in statis- 
tics for the deuteron and a wider coverage towards both 
lower and higher Q 2 and higher W . Once analyzed, this 
vastly larger data set will allow us to investigate in de- 
tail resonance electro-production on the neutron and the 
approach of the first moment of gf and g™ towards the 
GDH sum rule at the real photon point. 
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